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The statistics of isothermal hnes and loops of the Cosmic Microwave Background (CMB) 
radiation on the sky map is studied and the fractal structure is confirmed in the radiation 
temperature fluctuation. We estimate the fractal exponents, such as the fractal dimension Dc 
of the entire pattern of isothermal lines, the fractal dimension Dc of a single isothermal line, 
the exponent in Korcak's law for the size distribution of isothermal loops, the two kind of 
Hurst exponents. He for the profile of the CMB radiation temperature, and He for a single 
isothermal line. We also perform fractal analysis of two artificial sky maps simulated by a 
standard model in physical cosmology, the WMAP best-fit A Cold Dark Matter (ACDM) 
model, and by the Gaussian free model of rough surfaces. The temperature fluctuations of 
the real CMB radiation and in the simulation using the ACDM model are non-Gaussian, in 
the sense that the displacement of isothermal lines and loops has an antipersistent property 
indicated by He ~ 0.23 < 1/2. 

KEYWORDS: cosmic microwave background radiation, isothermal lines and loops, fractal analysis, self- 
affinity, Hurst exponent 



1. Introduction 

Our interesting feature of the Cosmic Microwave Back- 
ground (CMB) radiation discovered by Penzias and Wil- 
son^' is its uniformity over the full sky: the tempera- 
ture of radiation is uniform at a value higher than one 
part per thousand. This is now interpreted as evidence of 
the occurrence of the Big Bang and the rapid expansion 
of the universe during its very early stage. The recent 
development of measurements using artificial satellites 
has revealed, however, the existence of small but def- 
inite temperature fluctuations, from which cosmologists 
are expecting to determine the origin of large-scale struc- 
tures of galaxies and precise information for testing the 
Big Bang theory.^^ In the 1990s, the Cosmic Background 
Explorer (COBE) satellite conclusively detected a CMB 
radiation-temperature fluctuation in AT ~ 30/iK on an- 
gular scales larger than ~ 7°. 3. 4) PoUowing the launch of 
COBE, NASA constructed a second satellite, the Wilkin- 
son Microwave Anisotropy Probe (WMAP), which was 
launched on 30 June 2001. WMAP has made a map of 
the temperature fluctuation of the CMB radiation with 
much higher resolution and accuracy than C0BE.5) In 
the present paper, we discuss the temperature fluctuation 
of the CMB radiation from the viewpoint of statistical 
physics. In particular, after de Gouveia Dal Pino et al. 
analyzed the COBE data,^-* we represent structures of 
fluctuations observed by the WMAP by producing a set 
of isothermal lines and loops on the sky map. We extend 
the method of de Gouveia Dal Pino et al. by characteriz- 
ing fractal structures, and in the present paper not only 
the fractal dimension of single isothermal line but also 
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Fig. 1. Isothermal lines and loops of the CMB radiation on the 
sky map. Thick line represents one of the single isothermal lines. 



other fractal exponents of isothermal lines and loops are 
reported. 

2. Fractal analysis for isothermal lines and loops 

Figure 1 shows the set of isothermal lines and loops at 
an average temperature (2.725 K) in a partial region in 
the 5-year-CMB fluctuation measured by WMAP. The 
isothermal lines are produced from data of the WMAP 
Internal Linear Combination (ILG) Map provided at 
NASA's website. To characterize "spatial" fluctuations 
of the CMB radiation-temperature shown in Fig. 1, we 
evaluate the fractal dimensions for the entire set of 
isothermal lines and for a single isothermal line (a 
contour, for example, a thick line shown in Fig. 1) on the 
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CMB sky map. Moreover, we study the cumulative size 
distribution of a set of isothermal loops larger than x and 
show that it follows a power law with the characteristic 
exponent C,: 



n(x')dx' ~ X 



(1) 



This power-law distribution is known as Korcak's law.®' 

Now we report the results for the set of isothermal 
lines and loops at an average temperature. We analyze 
the 5-year-CMB radiation corresponding to the ILC and 
foreground-reduced W-band maps, covering ca. 25° on 
both sides of the line that exactly divides the galactic 
cast from the west in the 5-year-CMB sky maps measured 
by WMAP, taken horizontally.^' The foreground reduced 
map was removed using the foreground template method 
of Hinshaw et al.^°' and Page et al.-'^-'^) Moreover, we use 
rotated maps in which the north pole is centered to avoid 
the high contamination of the data at the Galactic plane. 
The fractal dimension D^. is determined using the box 
counting method. We scale the number of boxes N\^{e), 
covering the entire set of isothermal lines, as a function 
of the box size e as follows: 



^b(e) 



(2) 



Figure 2 shows the log-log plots of -/Vb(e) vs e on the ILC 
map. The slope in Fig. 2 yields the value of Dg as 



1.77. 



(3) 



On the other hand, the fractal dimension Dc is ob- 
tained by scaling the length s{Rg) of an isothermal 
line (a contour) as a function of the radius of gyration 



Ri 



(4) 



where is the center of mass of the isothermal line, 
is the position of the i-th pixel of an isothermal line and 
TVt is the total number of pixels included in the line, 
respectively. Figure 3 shows the log-log plots of s versus 
Rg on the ILC map. The slope in Fig. 3 yields Dc as 



1.39. 



(5) 



In the previous study, de Gouveia Dal Pino et al. ob- 
tained the fractal dimension Dc — 1.42 from the relation 
between the perimeter length and the area of an isother- 
mal loop generated from the CMB sky map measured by 
COBE-DMR.^' We also make a log-log plot of the cu- 
mulative size distribution of isothermal loops on the ILC 
map, as shown in Fig. 4. Figure 4 shows that Korcak's 
law (1) is established and 



C ~ 0.87. 



(6) 



The WMAP science working group reported that the ILC 
map provided a reliable estimate of the CMB signal on 
large angular scales greater than ^ 10°. i'^' On the other 
hand, on smaller scales, they also reported that the ILC 
map might not be as reliable as it should be. In order to 
check the validity of the fractal analysis for the ILC map, 
we reanalyzed the isothermal lines on the foreground- 
reduced W-band map. The obtained values are Dc — 



slope = 1.77 
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Fig. 2. Evaluation of Dc for the set of isothermal lines and loops 
on the CMB sky map. 
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Fig. 3. Log-log plots of the length s versus the radius of gyration 
Rg of a single isothermal line on the CMB sky map. 



1.74, Dc — 1.36 and C — 0.88. Because the differences are 
very small, we have concluded that the fractal analysis 
of the ILC map is reliable. 

Rough surfaces, such as landforms, are characterized 
by the Hurst exponent.®' ^'i^'^^-* The Hurst exponent H 
is defined by the scaling exponent of the PDF ph{s, h) 
of the fractional Brownian motion: 

1 { \ 



2s 



For iJo, we regard the "time" s as the horizontal position 
and the "displacement" h as the altitude. For He, we 
regard s as the length along the contour and h as the 
fluctuation of a contour. The characteristic exponents 
Dc , Dc , and C have scaling relations between them as 
reported by Matsushita et al.:^^' 



Dc 



Dc 



2 — He, 
2 

l + Hc 
Dc^2_ 
2 



(8) 
(9) 

(10) 
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Table I. Results of fractal analysis on the sets of isothermal lines. 
First column: for the real CMB sky map. Second column: for the 
artificial sky map simulated using the WMAP best-fit ACDM 
model. Third column: for the artificial sky map generated by the 
Gaussian white noise (Gaussian free model). Note that the mean 
values and standard deviation from 10 trials are shown for the 
artificial sky maps. 





CMB map 


ACDM model 


white noise 






1.77 


1.74 ± 0.03 


1.90 ± 0.02 






1.39 


1.36 ± 0.07 


1.53 ± 0.06 




c 


0.87 


0.8 ± 0.1 


0.9 ± 0.1 




from (8) 


0.23 


0.26 


0.10 




from (10) 


0.26 


0.4 


0.2 




from (11) 


0.22 


0.28 


-0.06 


He 


from (9) 


0.44 


0.47 


0.31 



Fig. 4. Cumulative area distribution of isothermal loops on the 
CMB sky map. 

where Hf, is the Hurst exponent of the profile of a given 
rough surface and He is the Hurst exponent of a single 
isoheiglit line or loop of the surface. In terms of the frac- 
tal dimension of a single isothermal line, Kondev and 
Henley derived the scaling relation^''^ 

Dc = (11) 

These scaling relations have been confirmed by a variety 
of data analyses of real landforms and numerical simula- 
tions of rough surface models. In general, when ffo = 1/2, 
the displacements at different points are independent of 
each other (as in usual Brownian motion), and when 
He > 1/2 (resp. H^ < 1/2) they have a positive (resp. 
negative) correlation, which is called persistence (resp. 
antipersistence)}^^^^^ The results for the CMB sky map 
are shown in the first column in Table I. The Hurst expo- 
nents for the CMB sky map evaluated by the four rela- 
tions eqs. (8)-(ll) are all less than 1/2. That is, the dis- 
placement of the isothermal lines and loops of the CMB 
radiation has antipersistent property, which is "noisier" 
than the usual Brownian motion. 

3. Numerical Simulations 

Chiang et al. performed a Monte Carlo simulation of 
the temperature fluctuations of the CMB sky map.^^-' 
They used the standard model in physical cosmology, 
which is called the WMAP best-fit A cold dark matter 
(ACDM) model. They have also provided an artificial sky 
map as reference for which the temperature fluctuation 
is generated from Gaussian white noise. It may corre- 
spond to the limit _ffo — > of the Gaussian rough sur- 
face model, which is called the Gaussian free fleld 
model in mathematical physics. For these two artifi- 
cial sky maps, we have also produced sets of isothermal 
lines and loops and performed the fractal analysis men- 
tioned above. The results are shown in the second and 
third columns in Table I. The obtained fractal dimen- 
sions Dc and for the ACDM model are very close to 
those for the real CMB sky map. In addition. Hurst expo- 
nents He and H^ obtained from the scaling relations eqs. 



(8)-(10) for the ACDM model is consistent with those 
of the CMB sky map. On the other hand, the artificial 
sky map generated from the Gaussian white noise has 
De ~2,De~ 1.5 = 3/2, C ^ 1, and H^ ~ 0. These val- 
ues arc expected in the Gaussian free field modcl^^'^"'^'^'' 
and are definitely different from the actual CMB values 
except C. 

4. CMB and statistical physics 

One of the well-discussed topics regarding the recent 
progress in the conformal field theories and statistical 
physics of critical phenomena and random fractal pat- 
terns is the introduction of the Stochastic Loewner Evo- 
lution (SLE) by Schramm.24) SLE is a stochastic differ- 
ential equation with a parameter, k, generating a random 
fractal curve on a plane. In particular, SLE with some 
characteristic values of k describes critical curves in sta- 
tistical physics such as the percolation exploration pro- 
cess and the phase boundary of the critical Ising model, 
as well as the continuum limits of the loop-erased ran- 
dom walk, the self-avoiding random walk (SAW), the 
random Peano curve, and other fractal curves. The 
set of isothermal lines generated from the Gaussian white 
noise has — 1.50 = 3/2 (the Gaussian free field 
model) as previously mentioned. This correponds to the 
SLE curve with n = 4.^'^) We note that the fractal ex- 
ponents Dc — 1.39 and He — 0.23, evaluated in the 
present study for the set of isothermal lines and loops 
on the CMB sky map, can be compared with the expo- 
nents of the phase boundary of the critical Ising model, 
Dc = 11/8 = 1.375 and He = 1/4. This corresponds to 
the SLE curve with k = 3.^^) Note that the SAW corre- 
sponding to the SLEg/s curve has the fractal exponents 
Dc = 4/3 = 1.3 and He = 1/3 0.3. 

In general, the Hurst exponent on the (s, /i)-plane is 
given by the self-affine exponents Vs and in each di- 
rection: ^'^^ 

H^^. (12) 

The exponents z/g and Vh are obtained using the asymp- 
totics s ~ N"^' and h ~ N'^'^, where s and h are the 
root-square-mean displacements of the s and h compo- 
nents of the position of each pixel in a curve, respec- 
tively, and N is the curve length between two arbitrary 
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points on the pattern. Formula (12) makes it possible to 
evaluate the Hurst exponents directly from the profile of 
the temperature fluctuation and each isothermal loop on 
the CMB sky map. Moreover, the measurement of the 
self-affine exponents will enable us to discuss the valid- 
ity of the scaling relations and the connection to other 
fractal interface models such as the Karder-Parisi-Zhang 
(KPZ) equation. Multifractal analysis is also useful 
in characterizing the distribution of nonuniform fiuctua- 
tions.'^^^ In a previous study, Diego et. al. have applied 
multifractal analysis to the study of the CMB tempera- 
ture fluctuation measured by COBE-DMR to calculate 
the generalized dimension Dq and f{a) spectrum. The 
CMB temperature fluctuation data measured by WMAP 
and PLANCK'^^^ satellites should be analysed using the 
multifractal method. The set of isothermal loops on the 
CMB sky map has the property showing that the size 
of a loop changes significantly owing to the change in 
the threshold temperature to draw isothermal loops. The 
relation between a threshold temperature and the size 
of a loop is similar to the relation between occupation 
probability and the size of a percolation cluster.^''') This 
similarity raises a "percolation problem on the CMB sky 
map" . These studies are currently under investigation by 
our collaborators. 

5. Summary 

We have investigated the fractal exponents of the 
entire set of isothermal lines, Dc of each single isother- 
mal line, and C, for Korcak's law of the CMB radiation- 
temperature fiuctuations measured by WMAP satellite. 
Moreover, we have determined the Hurst exponents 
and He that respectively characterize the self-affinity of 
the radiation-temperature profile and each isothermal 
line on the galactic coordinates obtainded by applying 
the scaling relations eqs. (8)-(ll). In the same way, we 
have also obtained the fractal exponents of two artificial 
sky maps generated using the WMAP best-fit ACDM 
model and the Gaussian free model of rough surfaces. 
The values of the various scaling exponents shown in Ta- 
ble I imply that 1) the fractal exponents D^, D^, ^e, 
and He of the real CMB sky map and those of the ACDM 
model coincide with each other, and 2) the evaluated He 
for the CMB radiation is less than 1/2, that is, the tem- 
perature fiuctuation has an antipcrsistent property. In 
cosmology, it is important to check whether the primor- 
dial fluctuation in the universe is Gaussian. To check 
for Gaussianity, the statistical properties of the CMB 
radiation-temperature fluctuation have been investigated 
by various methods, such as the use of angular bispec- 
^j.yj-i^35-39) Minkowski functionals,^^'^^"'^^-' phase map- 
ping techniques,"'^' and fractal analysis. ^^'^^' On the ba- 
sis of these analyses, it was found that the observations 
of the CMB radiation-temperature fluctuation are com- 
patible with the Gaussian hypothesis.^^, 35, 38, 39, 4i, 43) 

note, however, that Yadov and Wandelt presented evi- 
dence of the primordial non-Gaussianity /nl 0, where 
/nl denotes the nonlinearity parameter of the local 
type.3^) 

In the present study, we confirmed that the ensemble 
of isothermal lines and loops of the CMB has the fractal 



property and that the methods for analyzing random sur- 
faces and their isoheight lines and loops developed in the 
statistical mechanics and fractal physics are very useful. 
Here, we have reported that the fiuctuation of isothermal 
lines and loops of the CMB radiation has antipcrsistent 
property, since the Hurst exponent He < 1/2. We think 
that our present study will open a new field, where cos- 
mologists and statistical physicists can enjoy communi- 
cations with each other to establish better understanding 
of the universe. 
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